OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Characterization of the Altered Gene Expression Profile 
in Early Porcine Embryos Generated from 
Parthenogenesis and Somatic Cell Chromatin Transfer 

Chi Zhou\ John Dobrinsky^, Stephen Tsoi\ George R. Foxcroft\ Walter T. Dixon\ Paul Stothard\ 
John Verstegen^ Michael K. Dyck^* 

1 Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, Alberta, Canada, 2 International Center for Biotechnology, Minitube of 
America, Mount Horeb, Wisconsin, United States of America 



Abstract 

The in vitro production of early porcine embryos Is of particular scientific and econonnic Interest. In general, embryos 
produced from in vitro Assisted Reproductive Technologies (ART) manipulations, such as somatic cell chromatin transfer 
(CT) and parthenogenetic activation (PA), are less developmentally competent than in wVo-derived embryos. The 
mechanisms underlying the deficiencies of embryos generated from PA and CT have not been completely understood. To 
characterize the altered genes and gene networks in embryos generated from CT and PA, comparative transcriptomic 
analyses of in vivo (IVV) expanded blastocysts (XB), IW hatched blastocyst (HB), PA XB, PA HB, and CT HB were performed 
using a custom microarray platform enriched for genes expressed during early embryonic development. Differential 
expressions of 1492 and 103 genes were identified in PA and CT HB, respectively. In comparison with IW HB. The "elF2 
signalling", "mitochondrial dysfunction", "regulation of elF4 and p70S6K signalling", "protein ubiquitinatlon", and "mTOR 
signalling" pathways were down-regulated In PA HB. Dysregulatlon of notch signalling-associated genes were observed In 
both PA and CT HB. TP53 was predicted to be activated In both PA and CT HB, as 136 and 23 regulation targets of TP53 
showed significant differential expression In PA and CT HB, respectively, in comparison with IVV HB. In addition, 
dysregulations of several critical pluripotency, trophoblast development, and implantation-associated genes {NANOC, 
GATA2, KRT8, LGMN, and DPP4) were observed in PA HB during the blastocyst hatching process. The critical genes that were 
observed to be dysregulated in CT and PA embryos could be indicative of underlying developmental deficiencies of 
embryos produced from these technologies. 
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Introduction 

The domestic pig is not only an economically important 
livestock species, but is also an increasingly recognized biomedical 
animal model. Therefore, the in vitro production of early porcine 
embryos is of particular scientific and economic interest. Embryos 
produced from in vitro based systems using Assisted Reproductive 
Technologies (ART) are generally less developmentally competent 
in comparison with in vivo embryos. In swine, the in vitro production 
of pre-implantation embryos is much less efficient than in many 
other mammalian species (such as cattle) [1]. In vitro ART 
manipulations could have perturbing effects on embryonic gene 
expression, which potentially results in important negative long- 
term consequences [2], without displaying significant changes in 
tlie embryos' pre-implantation morphological characteristics [3- 
5]- 

Somatic cell nuclear transfer (SCNT) is a technology with great 
potential applications in basic and biomedical researches. How- 
ever, the application of SCNT is limited by its low embryonic 



survival rate and the high incidence of abnormalities in individuals 
that develop to term, and are believed to be associated with the 
incorrect or incomplete nuclear reprogramming [6,7]. Somatic 
cell chromatin transfer (CT) is a cloning technology that was 
designed to facilitate the reprogramming process [4,8], which 
involves in vitro remodelling of the donor nuclei prior to their 
transfer into enucleated oocytes to remove nuclear components 
that may interfere with nuclear remodelling [8]. Although 
promising results have been reported using chromatin transfer 
(CT), the CT-derived embryos still exhibit abnormalities similar to 
those observed following conventional SCNT [7,8]. Embryos 
derived from parthenogenetic activation (PA) are valuable for 
studies on gene imprinting [9] and are a potential alternative 
source of embryonic stem cells [9,10]. However, embryos 
generated from PA experience severe development failure [11]. 
The molecular mechanisms behind the deficiencies of embryos 
generated from PA and CT are not completely understood. 

The blastocyst is an embryonic stage that is frequently 
transferred into female recipients after ART manipulation 
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[12,13] in pig, and is, therefore, of particular scientific and 
economic interest. Similar to most mammalian species, the 
porcine blastocyst stage embryo has a distinct morphological 
structure that consists of the inner cell mass (ICM), internal cavity 
(blastocoele), and a single layer of epithelial trophectoderm (TE) 
with, or without (after hatching), the protective zona peUucida 
[14,15]. Blastocyst stage embryos need to hatch from the zona 
peUucida before implantation. Following blastocyst formation, the 
embryo expands in size and hatches from the zona pellucida to 
become a "free floating" hatched blastocyst in the uterus 
approximately 5—6 days after fertilization [16]. This "hatching" 
process is a critical and tightly regulated event during early 
embryonic development and any dysregulation of the hatching 
process leads to implantation failure and results in early embryonic 
loss [17]. 

In the present study, comparative transcriptomic analyses of in 
vivo (IW) expanded blastocysts (XB), IW hatched blastocyst 
(HB), PA XB, PA HB, and somatic ceU chromatin transfer (CT) 
HB were performed using a custom microarray platform enriched 
for genes expressed during early embryonic development - 
EmbryoGENE Porcine Array Version 1 (EMPVl, NCBI Gene 
Expression Omnibus (GEO): GPL14925) [18]. 

The objectives of the present study were (1) to characterize the 
effect of somatic cell chromatin transfer (CT) and parthenogenetic 
activation (PA) on the gene expression patterns of HB stage 
porcine embr^'os; (2) to identify critical genes and gene networks 
that were dysregulated during blastocyst hatching in PA embryos. 

Materials and Methods 

Animal ethics statement 

AU animal studies were conducted in accordance with the 
Canadian Council on Animal Care (CCAC) Guidelines and 
Policies with approval from the Animal Care and Use Committee: 
(Livestock) for the University of Alberta (Permit Number: DYCK- 
2006-56). 

Recovery of in vivo embryos 

In vivo (IW) derived porcine Germinal vesicle (GV), Mil, 2-ceU 
(2C). 4-ccII (4C), 8-cell (8C), morula(M), early blastocyst (EB), 
expanded blastocyst (XB), hatched blastocyst (HB) and embryonic 
day 11 (Dll) HB (hatched blastocyst before elongation) stage 
embryos were collected from gilts as described previously [19] and 
stored individually. The day of artificially insemination is 
considered day 0 (DO). All embryo samples were placed on dry 
ice immediately after collection and stored at — 80°C until RNA 
extraction. 

Production of in vitro-derived embryos 

All of the in vitro (somatic cell nuclear transfer (CT) and 
parthenogenetic activation (PA)) embryos used in the present study 
were produced by the International Center of Biotechnology, 
Minitube of America, MT Horeb, Wisconsin, USA (http://www. 
minitube.com). In brief, the CT reconstructed embryos were 
produced by using the Chromatin Transfer technology [20,21] 
under license from Hematech to Minitube (Verona, WI, USA). 
Oocyte' ('()lle[:tion, maturation, and micromanipulation were 
performed following established standard operating procedures 
[21,22]. The CT reconstructed embryos (for CT embryo 
production) and mature oocytes (for PA embryo production) were 
activated with incubation in 15 nM calcium ionomycin (Calbio- 
chem, CA, USA) supplemented mNCSU23 medium (Minitube, 
WI, USA) and subsequently an incubation of 1.9 mM 6- 
dimethylaminopurine (DMAP) supplemented mNCSU23 medium 



following previously established procedures [22]. The in vita 
activated CT reconstructed embryos and the parthenogenetically 
activated oocytes were both cultured in the PorcPRO mNCSU-23 

(Minitube, WI, USA) pig embryo culture medium system in 
38.7°C, 5% CO2, and 95-98% humidity for up to 8 days for 
expanded blastocyst and hatched blastocyst development. 

AU of the in vita (CT and PA)-derived embryo samples were 
placed on dry ice immediately after coUection and stored at 
-80°C untU RNA extraction. 

Total RNA isolation 

Total RNA was extracted from pools of 5 embryos using 
Arcturus PicoPure RNA Isolation Kit (Applied Biosystems, CA, 
USA). The RNA quality and integrity of each total RNA sample 
was evaluated by Bioanalyzer RNA 6000 Pico LabChip (Agilent 
Technologies, ON, Canada). Except for samples from Mil, 2C, 
4C, and 8C stage embryos, only high quality RNA samples (RNA 
integrity number (RIN) a7.5) were used for subsequent RNA 
amplification. It has been demonstrated that embryos from pre- 
embryo genome activation (pre-EGA) stages contains very low 
amounts of 28S rRNA which results in lower total RNA RIN 
value [23]. Therefore, total RNA samples from Mil, 2C, 4C, and 
8C stages with lower RIN values (range from 5.8 to 6.8, clear 188 
and 28S bands with no visual evidence of degradation) were 
utilized in this study. 

Microarray experimental design 

The comparative transcriptomic analyses were performed using 
a custom designed porcine embryo-specific microarray platform 
(EMPVl: EmiaryoGENE Porcine Array Versionl [GPL14925]) 
[18]. 

To characterize the effects of m vitro manipulations (PA and CT) 
on the porcine blastocyst transcriptome, comparative transcrip- 
tomic analyses among in vivo XB, in vivo HB, PA XB, PA HB, and 
CT HB were performed. Total RNA samples extracted from pools 
of 5 embryos from the same stage were amplified, labeUed with 
Cy5 dye, and hybridized with a Cy3 dye-labeUed reference 
amplified RNA (aRNA) pool on EMPVl microarray following a 
reference design [24-26] using three biological replicates from 
each group. 

Agilent two-colour RNA Spike-In (Agilent Technologies, ON, 
Canada) were amplified, labeUed and utilized as positive controls 
in each hybridization reaction as previously described [18]. 

Reference amplified RNA (aRNA) pool generation 

A reference aRNA pool was generated from 10 different 
embryonic stages (GV, Mil, 2C, 4C, 8C, M, EB, XB, HB, and 
Dll HB). Total RNA samples were amplified individually using 
RiboAmp HS^'"' kit (Apphed Biosystems, CA, USA); 1 ng of total 
RNA was used in each amplification and the quality and quantity 
of each aRNA sample was assessed using Bioanalyzer RNA 6000 
Nano LabChip (AgUent Technologies, ON, Canada) and Nano- 
drop ND-1000 (NanoDrop Technologies, Wilmington, USA). A 
total of 360 l^g of reference aRNA was generated by pooling 36 |ig 
amplified aRNA from each of the 10 embryonic stages. The 
reference aRNA pool was stored in ahquots at — 80°C until use. 

When applied to the EMPVl platform, the reference aRNA 
pool produced reference signals (signals that were higher than the 
average signal of negative controls) for 95% of aU the genes spotted 
on the microarray. 
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RNA amplification and labelling for microarray analysis 

Due to the low quantities of each total RNA sample, all RNA 
samples were amplified using RiboAmp HS^'"' kit (Applied 
Biosystems, CA, USA) following the manufacturer's instructions 
and generated amplified RNA (aRNA) targets for microarray 
reactions. One ng total RNA was utilized in each amplification 
reaction, and the quantity and quality of the aRNA products from 
RNA amplification reactions were evaluated by the Nanodrop 
ND-1000. Two |ig of aRNA were used in each labelling reaction. 
All labelling reactions were performed using the ULS Fluorescent 
Labelling Kit (Kreatech Diagnostics, Amsterdam, Netherlands) 
following the manufacturer's instructions. The labelling of aRNA 
targets was processed under an ozone-free environment. Probe 
concentration and labelling efficiency of each labelled sample was 
evaluated using Nanodrop ND-1000. 

Microarray hybridization, washing and data acquisition 

The hybridization, washing and drying steps of EMPVl 
microarray were conducted following the procedure described 
previously [18]. In short, aRNA samples were labelled with 
different dyes (Cy5 or Cy3) and hybridized on one microarray. 
Arrays were then incubated at 65 °C with rotation at 10 rpm for 
1 7 hours. After washing and drying steps that strictiy followed the 
instructions in Agilent manual, microarrays were immediately 
scanned using an Axon 4200AL scanner (Molecular Device, 
Sunnyvale, USA). The microarray results were submitted to NCBI 
Gene Expression Omnibus (GEO) Database (GSE48292). 

Microarray data were analys(;d using the FlexArray software 
package, which uses R and Bio-Conductor [27] and provides a 
user-friendly interface that facilitates data processing, visualization, 
and statistical analysis (Michal Blazejczyk, Mathieu Miron, Robert 
Nadon (2007). FlexArray: A statistical data analysis software for 
gene expression microarrays. Genome Quebec, Montreal, Cana- 
da, URL http://genomequebec.mcgill.ca/FlexArray). Simple 
background subtraction and within-array global loess normaliza- 
tion was performed on raw data from each array using the 
FlexArray software package. The threshold for positive spot 
selection from microarray data was determined as the mean value 
of all the negative control spots plus two standard deviations [18]. 
To identify differentially expressed genes, the normalized micro- 
array data was analyzed using the "limma" package [28] of Bio- 
conductor through FlexArray under the Benjamini and Hochberg 
false discovery rate (BH-FDR) [29] multiple comparison correc- 
tion condition through FlexArray [30]. For any particular 
comparison, only genes with a BH-FDR adjusted P value (B-H 
P-value) < 0.05 and a fold change (EC) > 2 (or <0.5) were 
considered to be significantiy up- or down-regulated. 

Gene expression data analysis 

Expression data obtained from the comparative transcriptomic 
analysis were analysed using the IPA (Ingenuity Pathway Analysis, 
Ingenuity Systems, www.ingenuity.com) Biological Functions 
Analysis, Canonical Pathway Analysis, and Upstream Regulator 
Analysis tools. The biological functions and canonical pathways 
analyses were performed under BH-FDR multiple testing correc- 
tion conditions. Only the biological functions and canonical 
pathways with a BH-FDR corrected P-value (B-H P-value) <0.05 
were considered significant. IPA Upstream Regulator Analysis 
predicts the activation status of the upstream regulator by 
calculating a regulation Z-score and an overlap P-value, which 
were based on the number of known regulation target genes from 
the dataset of interest, expression changes of these target genes, 
and their agreement with literature findings. Upstream regulators 
with an overlap P-value S 0.05 and an IPA activation Z-score S 



2.0 (or £ —2.0) were considered significantiy activated (or 
inhibited). Description of the calculation of the IPA regulation 
Z-score and overlap P-value is available in IPA white papers "A 
Novel Approach to Predicting Upstream Regulators". A full 
description of IPA analysis is available on the IPA website (http:/ / 
www.ingenuity.com) under "Upstream Regulator Analysis", 
"Biological Functions Analysis", and "Ingenuity Canonical 
Pathways Analysis". 

Real-time Quantitative PGR verification of gene 
expression results 

Fourteen genes selected from the comparative gene expression 
data were evaluated using SYBR Green I-based Real-time 

Quantitative PGR (QPCR). The primer sequences for all target 
genes are listed in Table 1 . A total of 1 ng total RNA isolated from 
each pool of 5 embryos was reverse transcribed into cDNA using a 
high capacity reverse transcriptase (Superscript VILO cDNA 
Synthesis Kit, Invitrogen) following the manufacturer's instruc- 
tions with a 20 |il reverse transcription (RT) reaction volume. An 
equal amount (10000 copies) of a synthetic RNA transcript (Xeno 
RNA Control, SYBR Green Cells-to-CT Control Kit, Ambion) 
was added to each reverse transcription reaction to serve as an 
external reference for SYBR Green I-based QPCR analysis, and 
as a positive control for reverse transcription, in order to assess 
variability resulting from any RT or PGR inhibitors. The cDNA 
products were then diluted 5 times, and 2.5 [il of the diluted 
cDNA was used as the template in each of the QPCR reactions 
performed with StepOnePlus Real-Time PGR System (Life 
technologies) and Fast SYBR Green Master Mix (Apphed 
Biosystems). The QPCR data was normalized with the external 
control gene (Xeno RNA Control, Ambion) using the qbase^"^^^ 
software (Biogazelle) [3 1] . The normalized QPCR data was then 
further analysed using the 2 ^^'-'^ method [31,32] to determine 
the relative cUfferential expression (fold changes) of each target 
gene. 

Results 

Altered gene expression profile in PA- and CT-derived HB 

The reference design, which was used in the microarray 
comparative transcriptomic analysis among embryos derived from 
the PA and CT, allows for reliable comparisons among different 
groups in the analysis [26], as described in the methods. 

Comparative microarray analysis revealed 1492 and 103 
significant differentially expressed (FC > 2 or < 0.5, B-H P- 
value<0.05) genes in PA- and GT-derived HB, respectively, in 
comparison with IW-derived HB (Dataset SI). In comparison 
with IW HB, 55 genes showed significant differential expression 
in both PA and CT HB, and 54 out of these 55 genes showed the 
same direction of expression changes (up- or down-regulation) in 
PA and CT HB. 

IPA biological function (bio-function) analysis revealed 
19 and 48 biological function categories that were significantiy 
altered (B-H P-value<0.05, and have more than 8 molecules 
included in the analysis) in PA- and CT-derived HB, respectively 
(Dataset S2). The most significantiy altered (B-H P-value<0.01, 
and have mor(; than 8 molecules included in the analysis) bio- 
function categories in PA HB and CT HB were further identified 
(Figure 1). The four most significantiy altered bio-functions in PA 
HB were associated with "cellular growth and proliferation", 
"cellular development", "cell cycle", and "neurological disease"; 
and the four most significantly altered bio-functions in CT HB 
were associated with "cell cycle", "neurological disease", "skeletal 
and muscular disorders", and "nucleic acid metabolism". The 
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Table 1. Primer sequences used in Real-time PCR verification. 




Gene symbol 


Associsted Porcine 
RefSeq Accession No. 


Primer 


Primer sequence (5'-3') 


Product size 
(bp) 


LGMN 


XM_001 927082 


Forward 


AGACGCTCCACAAACAGTAC 


95 






Reverse 


CAACTTCATGGCAGAGATGGA 




GATA2 


NM_213879 


Forward 


CTCCAGCTTCACCCCTAAG 


157 






Reverse 


CCCGTTCATCTTGTGGTACAG 




KRT8 


NM_001159615 


Forward 


AGATCCAAAAGCGTACCGAC 


136 






Reverse 


AGCTGCCTGTAGAAGTTGATC 




PSEN2 


NM_001 078666 


Forward 


CTCAACTCCGTGCTCAACA 


148 






Reverse 


GATGTAGGTGAAGAGGAAGAGC 




NCSTN 


XM_001 928786 


Forward 


CCCCGCAATGTCATGTTTG 


92 






Reverse 


AACTTGCCCTTCTCCATATCG 




HESl 


NM_001 195231 


Forward 


CTGGAGAAGGCGGACATTC 


92 






Reverse 


GCTCGGGTCTGTGCTTAG 




HEY2 


NM_001 243329 


Forward 


CTGCAAAGTTAGAAAAGGCCG 


145 






Reverse 


TCTGTTAAGCACTCTCGGAATC 




ANXA8 


NM_001 243599 


Forward 


AGACATACAAGCAGATACCAGTG 


142 






Reverse 


CTTCTCACCCGCTGCATAC 




SLC36A2 


XM_0031 34141 


Forward 


CATCACCCAGTACATCATCCAG 


127 






Reverse 


CAGAACCACACCAATGCTTTC 




KCTD3 


XM_003357619 


Forward 


AGAAGTTCCCTCTGCGAATG 


149 






Reverse 


CGTACCATAGGCGATCTCAATC 




NANOG 


NM_001 129971 


Forward 


GGACTTTTCCTACAATCCAGC 


153 






Reverse 


CCCATAAACCTCAGGCATTG 




JAGl 


XM_001 926559 


Forward 


ACATAGCCCGAAACAGTAGC 


158 






Reverse 


GTTGTAGCAGGGATGAGGAC 




DPP4 


NM_214257 


Forward 


TGCGGATTCCATACCCAAAG 


137 






Reverse 


ATCCCCTATTAACACAGACGC 




KRT18 


XM_003126180 


Forward 


TTGACCGTGGAGTTGGATG 


149 






Reverse 


ACCACTGAGGTGCTCTCC 




Xeno 


Control primer Xeno from SYBR Green Cells-to-CT Control Kit (Ambion) 


105 


doi:l 0.1 371 /journal.pone.0091 728.t001 



"cell cycle", and "neurological disease"-associated bio-functions 
were significantly altered in both PA and CT-derived HB stages 
embryos. 

IPA canonical pathway analysis revealed eight canonical 
pathways that were significantly altered (B-H P-value<0.05, and 
have more than six molecules included in the analysis) in PA HB in 
comparison with IW HB (Figure 2). The "eIF2 signalling", 
"mitochondrial dysfunction", "regulation of eIF4 and p70S6K 
signalling", "protein ubiquitination", and "mTOR signalling" 
pathways were the five most significantly changed canonical 
pathways between PA HB and IW HB. Specifically, most of the 
differentially expressed genes associated with these pathways were 
down-regulated in PA HB (Figure 2). 

IPA upstream regulator analysis revealed five (MYC, 
MYCN, NOBOX, PPARGCIA, and TP53) and one (TP53) 
transcription factors that predicted to be significantly activated (or 
inhibited) in PA and CT-derived HB, respectively, in comparison 
with IW HB (Dataset S3). Transcription factors PPARGCIA, 
MYC, and MYCN were predicted to be inhibited in PA HB; and 
transcription factors NOBOX and TP53 were predicted to be 
activated in PA HB. 



Transcription factor TPS 3 was predicted to be significandy 
activated in both PA and CT HB. Although no significant 
differential expression of the TP53 gene was observed, 136 and 23 
regulation targets of TPS 3 showed significant differential expres- 
sion in PA and CT-derived HB, respectively, in comparison with 
IW embryos (Dataset S4). In addition, 1 1 regulation targets 
{ANXA8, CTSH, CTSK, GSTPl, HSP90AA1, IL6, Mr06, PERP, 
PHLDA3, PRDX3, and PSEN2) of TP53 showed differential 
expression in both PA and CT HB compared with IW HB. 
The down-regulation of PSEN2 (Figure 3A) and the up-regulation 
of ANXAS (Figure 3B) in PA and CT HB were confirmed by 
QPCR. ANXA8 displayed detectable expression levels in both PA 
and CT HB, and AMXA8 displayed significantly higher expression 
in PA HB than CT HB. No detectable expression oi ANXA8 was 
observed in IW HB by QPCR analysis. 

Significant differential expression of four "notch signaUing"- 
associated genes {PSm2, HEr2, HESl, and JAGl) were observed 
in PA HB in comparison with IW-derived HB embryos. In 
comparison with IW HB, the microarray analysis revealed 
significant down-regulation ot HET2, HESl, and JAGl genes, and 
significant up-regulation of PSEN2 showed in PA HB. Significant 
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IPA Bio-function Category 



Figure 1. Altered biological function categories in PA and CT-derived HB. Bar chart shows the significantly altered biological function 
categories in IPA biological function (bio-function) analysis. Major Y axis on the left shows the number of differentially expressed genes that involved 
in the biological function category. Secondary Y axis on the right shows the significance (-log (B-H P-value)) of the altered biological function 
category. The orange line shows the significance threshold of cut off of -log (B-H P-value = 0.05). 
doi:1 0.1 371/journal.pone.0091 728.g001 
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IPA Canonical Pathways 



■ number of genes down-regulated in PA HB 
(PA HB vs. IW HB) 

■ number of genes up-regulated in PA HB (PA 
HBvs. IWHB) 

■ number of genes down-regulated in CTHB 
(CTHBvs. IWHB) 

■ number of genes up-regulated in CT HB (CT 
HBvs. IWHB) 

- -log(B-H p-value) PA HB vs. IW HB 

- -log(B-H p-value) CT HB vs. IW HB 

- Threshold -log(B-H p-value=0.05) 



Figure 2. Altered canonical pathways in PA and CT-derived HB. Bar chart shows the altered canonical pathways in IPA canonical pathways 
analysis. Major Y axis on the left shows the number of differentially expressed genes that involved in the canonical pathway. Secondary Y axis on the 
right shows the significance (-log (B-H P-value)) of the canonical pathway. The orange line shows the significance threshold cut off of -log (B-H P- 
value = 0.05). 

doi:1 0.1 371 /journal.pone.0091 728.g002 
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Figure 3. QPCR verification result. QPCR verification result of 14 selected genes. The mRNA expression levels of these genes were normalized 
with the external control gene (Xeno), and were calculated with 2^'^'^'"' relative quantification. Bar charts showing the relative expression levels of 
PSEN2, ANXA8, HES1, JAG1, HEY2, NCSTN, KRT18, KRT8, GATA2, NANOG, SLC36A2, KCTD3, DPP4, and LGMN genes in IVV XB, IW HB, PA XB, PA HB, 
and CT HB (KCTD3, SLC36A2, and LGMN genes were not tested in CT HB). The relative expression levels of in each sample were standardized with 
their expression Error bars shows the standard error (*: P < 0.05). Dashed lines indicate 1.0 expression level. ND: not detected. NT: not tested. 
doi:1 0.1 371/journal.pone.0091 728.g003 



down-regulation of the PSEN2 genes was also observed in CT- 
derived HB in comparison with IW HB. Another three [NCSTN, 
HESl, and JAGl) "notch signaUing"-associated genes showed 
altered expression in CT HB in comparison with IW HB, but 
with less statistical significance (FC>2 or <0.5, P-value<0.0.5 but 
B-H P-value >0.05). 

Five "notch signalluig"-associated genes (PSEN2, HEY2, HESl, 
NCSTN, and JAGl) were selected for QPCR verification. HESl 
(Figure 3C) was up-regulated in both PA and CT HB in 
comparison with IW HB. PSEN2 (Figure 3A) was down-regulated 
in CT HB in comparison with IW HB, and did not display 
detectable expression in PA HB. JAGl (Figwe 3D) did not display 
detectable expression in IW HB, but was expressed in both PA 
HB and CT HB. HET2 expression (Figure 3E) was not detectable 
in PA, CT, and IW HB embryos in the QPCR analysis. No 
significant expression change of NCSTN (Figure 3F) was observed 
among PA, CT and IW-derived HB in the QPCR analysis. 

In comparison with IW HB, significant down-regulation 
(FC = 0.3, B-H P-value<0.05) KRT18 (Figure 3G) was observed 
in PA HB, and a less significant down-regulation (FC = 0.66, P- 
value<0.05 but B-H P-value>0.05) oi KRTl 8 was observed in CT 
HB. QPCR analysis oiKRTlS (Figure 3G) and JCRTS (Figure 3H) 
expression showed that the KRT18 and KRT8 genes were down- 
regulated in PA HB, and the KRTIS was down-regulated in CT 
HB, in comparison with IW HB. 

In addition, microarray analysis revealed significant down- 
regulation (FC<0.5, B-H P-value<0.05) oiGATA2 and jWOG in 
PA HB in comparison with IW HB. QPCR analysis results 
confirmed this down-regulation of GATA2 in NANOG in PA HB 
(Figure 31 and 3J). However, the down-regulation o{ GATA2 in CT 
HB was not statistically significant. In the QPCR analysis, NANOG 



expression was only detectable in IW HB, and no detectable 
expression of NANOG was observed in PA and CT HB. 

Microarray analysis also revealed significant up-regulation (FC 
> 2, B-H P-value<0.05) of four precursor-microRNAs (pre- 
miRNA) {MIR1343, MIR149, MIR505, and MIR192) in PA HB in 
comparison with IW HB. Only trends (P-value<0.05 but B-H P- 
value>0.05) of differential expression of the pre-miRNA of 
MIR505 (FC = 0.57) and MIR192 (FC= 1.57) were observed in 
the CT HB in comparison with IW HB. 

Altered gene expression-regulation during blastocyst 
hatching of PA-derived embryos 

Comparative transcriptomic analysis among IW XB, IW HB, 
PA XB, and PA HB revealed that during the transition from XB to 
HB, differential expression (FC > 2 or <0.5, B-H P-value<0.05) 
of 3 and 3 1 genes were observed in PA and IW-derived embryos, 
respectively (Dataset S5). 

The comparative microarray analysis revealed three genes 
[KCTD3, ANNAS, and SLC36A2) that showed statisticaUy signifi- 
cant up-regulation from XB to HB in PA embryos. However, no 
significant differential expression of these three genes was observed 
in between IW-derived XB and IW HB. 

QPCR analysis confu-med the up-regulation of SLC36A2 
(Figure 3K) and ANXA8 (Figure 3B) from the XB to HB stage in 
PA embryos. SLC36A2 showed no significant differential expres- 
sion between IW XB and IW HB, and ANXA8 expression was 
not detectable in IW XB and IW HB. No significant differential 
expression of KCTD3 was observed between PA XB and PA HB in 
the QPCR analysis (Figure 3L). 

Significant up-regulation (FC > 2, B-H P-value<0.05) oi DPP4 
and LGMN from XB to HB in in aroo-derived embryos were 
observed in the microarray analysis. Trends toward up-regulation 
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of the DPP4 and LGMN were also observed in the PA embryos 
from XB to HB. In addition, a trend (P-value<0.05 but B-H P- 
value >0.05) of up-regulation of the trophectoderm development- 
associated gene KRT8 (FC = 1.9) and GATA2 (FC = 2.4) from XB 
to HB in IW embryo was observed in the microarray analysis. 
However, no differential expression of KRT8 was observed 
between PA XB and PA HB embryos. Results from QPCR 
analysis confirmed the up-regulation oiDPP4, LGMN, GATA2 and 
KRT8 from XB to HB in vim (Figure 3M, 3N, 3H, and 3G). In 
comparison with IW embryos, DPP4 and LGMN displayed a 
smaller up-regulation from XB to HB in PA embryos. No 
differential expression o{GATA2 and KRTS was observed between 
PA XB and PA HB in the QPCR analysis. 

Three (HEY2, HESl, and JAGl) "Notch signaUing"-associated 
genes showed down-regulation (FC<0.5), but with reduced 
statistical significance (P-value<0.05 but B-H P-value >0.05), 
from XB to HB in IW embryos in the microarray analysis. HESl 
showed more than 2.5 fold down-regulation from XB to HB in 
both IW and PA embryos. HET2 and JAGl showed more than 
2.4 fold down-regulation from XB to HB in IW embryos, but no 
significant differential expression of these two genes was observed 
in PA embryos. 

Results from QPCR analysis confirmed the up-regulation of 
HESl and the down-regulation of HET2 and JAGl from XB to 
HB in IW embryos (Figure 3C-E). Although up-regulation of 
HESl and down-regulation of JAGl from XB to HB in PA 
embryos were observed in the QPCR analysis, the expression 
changes of these two genes were less significant than 1 W embryos 
(Figure 3C-D). HEY2 displayed a higher expression in IW XB 
than PA XB, and HET2 expression was not detectable in both PA 
and rW HB in the QPCR analysis (Figure 3E). 

Discussion 

The embryos generated after in vitro manipulations such as 
parthenogenetic activation and nuclear transfer displayed slower 
and less effective development [33-36], and dysregulation of 
critical gene networks is probably associated with these deficien- 
cies. 

The first objective of the present study was to characterize the 
effects of somatic cell chromatin transfer (CT) and parthenogenetic 
activation (PA) on the gene expression patterns of hatched 
blastocyst stage porcine embryos. 

Comparative microarray analysis revealed 1492 and 103 
significantiy differentially expressed genes in PA and CT-derived 
HB, respectively, in comparison with IW-derived HB. This large 
gene expression profile differences between PA HB and IW HB 
observed in the present study is consistent with previous studies in 
different species [9,33-35]. The gene expression profile differences 
between CT and FW-derived HB observed in the present study 
was less pronounced than the differences previously reported 
between SCNT and IW-derive porcine blastocyst stage embryos 
[36]. 

In comparison with IW HB, the"eIF2 signalling", "mTOR 
signalling", "regulation of eIF4 and p70S6K signaUing", "mito- 
chondrial dysfunction", and "protein ubiquitination pathway" 

pathways were the 5 most significantly altered pathways in PA HB, 
and most of the differentially expressed genes associated with these 
5 pathways were down-regulated in PA HB. 

Eukaryotic translation initiation factor 2 (eIF2) plays a key role 
in the recognition of the correct start codon during translation 
initiation process [37]. Phosphor\lation of eIF2 reduces global 
translation and activates the transcription of "stress recovery" 
genes in response to environmental stresses such as amino acid 



deficiency, hea\ \ metal toxicity, and bacterial infection [37,38]. It 
has been reported that cells with defective eIF2 signalling were 
more susceptible to bacterial invasion [38]. The "mTOR 
signalling pathway" plays a critical role in the regulating of cell 
growth, proliferation, translation, protein synthesis and survival 
[39-41]. The eIF4 initiation factors are responsible for recruiting 
mRNA to a ribosome during translation process [42]. The 
translation eIF4 initiation factors and p70 S6 kinase (p70S6k) both 
play critical roles in the translation and protein synthesis 
regulation, and both eIF4 and p70S6k are regulation targets of 
mTOR [41,42]. Many emlronmental stimuli induding growth 
factors, hormones, and nutrient availability can regulate the eIF4 
and p70S6K through "mTOR signalling pathway" [39]. The 
down-regulation of genes associated with the "eIF2 signalling", 
"mTOR signalling", "Regulation of eIF4 and p70S6K signalling" 
pathways suggest that the general translation and protein synthesis 
are affected in PA HB; and many "mTOR signaUing"-associated 
critical biological processes are also significantiy affected in PA 
HB. 

Mitochondria, especially as an ATP generation source, are 
critical for the development of early embryos, and perturbation in 
their functions is associated with compromised embryonic 
competence [43] . Mitochondrial dysfunction in oocytes is directiy 
responsible for the high levels of developmental retardation and 
early arrest of pre-implantation embryos produced in vitro [44]. In 
the present study, the down-regulation of "mitochondrial 
dysfunction"-associated genes in PA HB suggests compromised 
mitochondria function in PA IIB. 

The "Ubiquitin-proteasome pathway" is responsible for the 
selective degradation of soluble cellular proteins in most cases [45]. 
Ubiquitination of cellular protein is essential for the ubiquitin- 
proteasome pathway-dependent cellular protein degradation [46] . 
Degradation of maternal proteins through the "ubiquitin-protea- 
some pathway" is believed to be important for the oocyte- to- 
embryo transition [47]. In this study, signiiicant differential 
expressions in genes associated with "protein ubiquitination 
pathway" were observed, suggesting an altered protein degrada- 
tion process in PA embr\'OS. 

TP53 (tumor protein p53) encoding a well- known cell-cycle 
regulator and apoptosis mediator [48], and it has been previously 
reported that the embryos derived from parthenogenetic activa- 
tion experience a higher apoptotic cell death rate [1 1] . Results 
from the present study showed that the TPS 3 is predicted to be 
activated in both PA and CT HB in comparison with the IW HB, 
where the number of differentially expressed TP53 regulation 
targets in PA HB was more than four times higher than the 
number of differentially expressed TP53 regulation targets in CT 
HB. In addition, AMXA8 (annexin A8) is a member of the annexins 
(ANXs) family, which is a group of Ca^^-dependent phosphohpid- 
binding proteins. ANXs are involved in many important biological 
processes including vesicle trafficking, calcium signalling, cell 
growth, cell cycle, and apoptosis [49]. 0\cr expression oi A]SXA8 
has been reported to be associated with cancer and apoptosis [50] . 
In the present study, ANXA8 displayed significantiy higher 
expression in PA HB than CT HB, and no detectable expression 
ofANXAS was observed in IW HB. These results suggest that an 
activated apoptotic process might be induced in both PA and CT 
derived HB, and that the activation of this apoptotic process 
appears to be greater in PA HB than in CT HB. 

NOTCH is an important regulator of development in many 
animals [51], which participate in many critical biological 
processes including cell fate specification, differentiation, prolifer- 
ation, apoptosis, migration, and angiogenesis [52]. Small pertur- 
bations in Notch activity could lead to numerous developmental 



PLOS ONE I www.plosone.org 



7 



March 2014 | Volume 9 | Issue 3 | e91728 



Embryo Transcriptome and Reproductive Technologies 



defects and diseases [51]. Notch signalling is initiated through 
ligand-receptor interactions between neighbouring cells [52]. The 
NOTCH-mediated HESl expression plays an important role in 
the regulation of cell fate decision [52]. In mammals, the two 
highly homologous presenUin genes [PSENl and PSEM2) play- 
important roles during early embryonic development and both of 
the presenilin genes are positive regulators of the "notch 
signalling" pathway [53,54]. Results from the present study 
showed that one of the mammalian Notch ligands Jagged 1 
(encoded by JAGl gene) [52,55] and two other members [HESl 
[56] and PSKN'2 [54,57]) of "Noti:h signalling" pathway were 
significantly differentially expressed between PA and IW-derived 
HB. Less dramatic differential expression of these three "notch 
signalIing"-associated genes were also observed in the CT HB. 
These results suggest that the "notch signalling" pathway is 
dysregulated in both PA and CT HB, and this dysregulation is 
more significant in PA HB than in CT HB. The altered regulation 
in Notch signalling probably contributes to the impaired 
development of PA and CT-derived embryos. 

As one of the key regulators of pluripotency, the transcription 
factor MANOG functions as a repressor of the extra-embryonic 
endoderm (ExE) or primitive endoderm (PE) cell fate [58]. In 
comparison with IVV HB, significant down-regulation oi NANOG 
in both PA and CT HB was observed in the present study, which 
suggests a compromised regulation of cell fate specification and 
TE differentiation in PA HB and CT HB. 

Transcription factor GATA binding protein 2 {GATA2) is 
expressed in trophoblast giant cells and acts as important regulator 
for trophoblast-specific gene expression and placental function 
[59,60]. Expression of GATA2 genes is essential for normal 
embryonic development in rodents [59]. Keratins 8 (KRT8), 
keratin 18 [KRT18) and keratins 19 (KRT19) are predominandy 
expressed in epithelial components of glandular tissues in rodents 
and humans [61—63]. Expression of keratin 8 and keratin 18/19 
are expressed in TE and are essential for the integrity of a 
specialized embryonic epitheUum (trophoblast giant cells) layer 
and the survival of embryos [60,62,63]. In the present study, 
GATA2, KRT18, and KRT18 showed significant down-regulation 
in PA HB, but only KRT18 showed significant down-regulation in 
CT HB embryos in the QPCR analysis. These results suggest 
impaired trophoblast development in both PA HB and CT HB, 
and trophoblast development in CT HB is less affected than PA 
HB. 

Although the DPP4 (dipeptidyl peptidase 4) was reported to be 
differentially regulated in the CT-derived bovine day 45 placenta 
[7] , no significant differential expression of DPP4 was observed in 
PA and CT HB in the present study. 

MicroRNAs (miRNA) are believed to be key regulators in pre- 
implantation embryonic development and difierentiation [64,65]. 
Recent reports suggest that the microRNA reprogramming is 
inc:omplete and inconsistent in cloned embr^'os [35,66]. In the 
present study, microarray analysis revealed significant differential 
expression of four pre-miRNAs in PA HB in comparison with IW 
HB. Two of these 4 pre-miRNAs showed trends of differential 
expression, and no statistically significant differentially expressed 
pre-miRNA was observed between CT and IW HB. During pre- 
implantation development of embryos, dynamic synthesis and 
degradation of miRNAs coexists [65]. Hence the differential 
expression of pre-miRNA does not guarantee the dififerential 
expression of mature miRNA. 

In this study, the oocyte in vitro maturation, in vitro activation, as 
well as the embryo in vitro culture processes for the CT and PA 
embryo generation followed exactly the same procedures. The 54 
common differentially expressed genes that were observed in CT 



and PA embryos in comparison with IW embryos could be 
associated with any of these common in vitro manipulation 
processes. Further studies are necessary before the differential 

expression of these "common differentially expressed genes" could 
be connected with any specific in vitro process. 

Tlu- second objective of the present study was to identify 
dysregulated genes and gene networks in PA embryos during 
blastocyst hatching. Hatching is a critical and necessary process 
during the early development of mammalian embryos. Blastocyst 
hatching is a well programmed and tightiy regulated event, and 
dysregulation of this critical process leads to implantation failure 
and results in early embryonic loss [17]. Dysregulation of critical 
genes and gene networks during blastocyst hatching process are 
probably contributed to the deficiencies in embryos generated 
from PA. 

In the present study, significant differential expression of 31 
genes were observed during the blastocyst hatching process (from 
XB to HB) in IW embryos, but these 3 1 genes were not properly 
regulated in PA embryos during blastocyst hatching. On the other 
hand, SLC36A2 and AJ^XAS showed significant up-regulation 
during the blastocyst hatching process in PA embr)'OS, but no up- 
regulation of these two genes were observed in IW embryos. 
SLC36A2 (Solute carrier family 36 (proton/ amino acid sympor- 
ter), member 2) mediates the transport of amino and fatty acids, 
which are critical to early embryonic development [67,68]. 
Further work is necessary to determine if this up-regulation of 
SLC36A2 is compensating for the function of other dysregulated 
genes in PA embryos and reflecting the increased need for 
nutrients in the rapidly developing embryos. 

LGMN (legumain), also known as cysteine protease 1, is 
involved in protein processing and is highly expressed in the 
placenta of pig [69]. Legumain has been reported to be expressed 
in bovine trophoblast and associated with the regulation of 
trophoblast invasiveness and endometrial remodelling during 
implantation [70]. DPP4 (dipeptidyl peptidase 4) is a membrane- 
bound aminopeptidase, which is associated with placental 
development and the establishment of proper fetal-maternal 
interactions in cattle and human [7,71]. In the present study, 
dramatic up-regulation of LGMN and DPP4 were observed during 
hatching process in IVV embryos, but the expression changes of 
Z,G/WV and DPP4 observed during hatching process in PA embryos 
are much less dramatic. Results from the present study showed 
that the expression of sc\-(^ral critical pluripotency, trophoblast 
development, and implantation-associated genes [NANOG, GATA2, 
KRT8, LGMN, and DPP4) were not properly regulated during the 
blastocyst hatching process in PA embryos. In addition, altered 
regulation of "notch signaIling"-associated genes was also 
observed during the blastocyst hatching process in PA embryos. 
Faihng to regulate the expression of these critical genes during the 
hatching process is probably contributed to the delayed and less 
efficient development of PA embryos. 

Further protein expression level data for these dysregulated 
genes that were identified in the present study would be necessary 
before a definite link could be drawn between dysregulations of 
these genes and deficiencies that observed in PA and CT embryos. 

Conclusion 

In the present study, we have successfully characterized the 
altered gene expression profiles in porcine HB embryos derived 
from parthenogenetic activation and somatic cell chromatin 
transfer, in comparison with in OT»o-derived HB. Specifically, we 
have identified several signalling pathways, critical genes, and 
critical gene networks that were significantiy altered in the PA- and 



PLOS ONE I www.plosone.org 



8 



March 2014 | Volume 9 | Issue 3 | e91728 



Embryo Transcriptome and Reproductive Technologies 



CT-derived HB stage embryos. In addition, we have also 
identified several critical genes that were not properly regulated 
during the blastocyst hatching process in embryos derived from 
PA. 

To date, morphological characteristics and blastocyst formation 
rat(" are still two of the major parameters commonly used in 
embryonic developmental competence assessment [2]. Results 
from the present study showed that embryos produced from PA 
and CT could develop into expanded blastocyst and hatched 
blastocyst stage, even with dysregulations of critical pathways and 
gene networks. Hence, the morphological criteria and blastocyst 
development ratio are insufficient to determine the ultimate 
competence of embryos generated after in vitro ART manipulations 
(such as PA and CT). The critical genes that exhibited altered 
expression in CT and PA embryos could be indicative of 
underlying developmental deficiencies of embryos produced from 
these technologies. 

Supporting Information 

Dataset SI Expression data of significant difierentially 
expressed (B-H P-value<0.05, F02 or <0.5) genes in 
PA HB vs. IW HB and CT HB vs. IW HB analyses. 

(XLSX) 

Dataset S2 IPA bio-fiinction analysis result of signifi- 
cantly altered (B-H P-value<0.03, molecules involved in 
the analysis £ 8) biological function categories in PA HB 
vs. IW HB and CT HB vs. IW HB analyses. 
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